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Введение. При математическом конечноэлементном моде- 
лировании используется усредненное значение механиче- 
ских характеристик материала деформируемых твердых 
тел. В авиа-, машиностроении, строительстве, медицине и 
других областях все шире применяются полимерные ком- 
позитные материалы и материалы природного происхож- 
дения. В последнем случае реальное изменение механиче- 
ских характеристик значительно отличается от усреднен- 
ного, следовательно, при использовании усредненных 
параметров для построения и анализа конечноэлементных 
моделей результаты могут существенно искажаться. 

В данной статье описано создание математических мето- 
дов исследования изменения механических характеристик 
материала неоднородных деформируемых твердых тел. 
Полученные таким образом результаты применены для 
построения конечноэлементных моделей и анализа их 
напряженно-деформированного состояния. 

Материалы и методы. В качестве неоднородных дефор- 
мируемых твердых тел рассмотрены материалы природ- 
ного происхождения и композиты. Для исследования из- 
менения механических характеристик материала разрабо- 
тан способ, основанный на использовании двух составля- 
ющих: пиксельной характеристики растровых изображе- 
ний сканирования компьютерным томографом и экспери- 
ментальных данных натурных испытаний стандартных 
образцов. 

Результаты исследования. Создан комплекс математиче- 
ских методов моделирования интерпретации растровых 
изображений сканирования компьютерным томографом, 
позволяющий проводить исследование любых сложных 
структур реальных деформируемых твердых тел. Резуль- 
таты используются при построении конечноэлементных 
моделей таких тел с учетом неоднородности механиче- 
ских характеристик материала. 

Анализ напряженно-деформированного состояния конеч- 
ноэлементных моделей тестовых образцов доказал точ- 
ность и сходимость численного решения метода конечных 
элементов при моделировании свойства неоднородности 
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механических характеристик материала. 

Обсуждение и заключения. Разработанный подход может 
быть применен для любых физических принципов скани- 
рования (рентгеновский, ультразвуковой, лазерный и др.) 
и для любых типов материалов, если информация, полу- 
ченная в результате сканирования, сформирована в виде 
цифрового (растрового) изображения. 


Ключевые слова: метод конечных элементов, деформи- 
руемое твердое тело, неоднородность, механические ха- 
рактеристики материала. 
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Введение. При конечноэлементном (КЭ) моделировании механические характеристики материала де- 
формируемых твердых тел (ДТТ) задаются в виде усредненного значения, используемого для всей модели. 
Например, механические характеристики стали определяются при испытании на растяжение и сжатие стан- 
дартных образцов правильной формы. Такой подход приемлем для расчетно-инженерной практики. Однако в 
современном авиа-, машиностроении, строительстве, медицине и других сферах все шире применяются поли- 
мерные композиты и материалы природного происхождения. Реальное изменение механических характеристик 
таких материалов сильно отличается от усредненного значения, следовательно, при использовании усреднен- 
ных параметров для построения и анализа конечноэлементных моделей результаты могут существенно иска- 
жаться. 

Результаты исследования 

1. Основная математическая зависимость моделирования неоднородности реальных деформиру- 
емых твердых тел. Для решения представленной проблемы в работе предлагается способ определения реаль- 
ного изменения механических характеристик материала, основанный на использовании: 

— результатов сканирования ДТТ; 
— данных натурных испытаний стандартных образцов [1]. 

Такой способ дает возможность распознавать любую сложную неоднородную структуру материала и 
применить эти данные для повышения точности и реалистичности при математическом конечноэлементном 
моделировании реальных ДТТ. 

На рис. 1 а представлен пример растрового изображения сканирования материала природного проис- 
хождения (человеческой кости). В ее сечении видно изменение механических характеристик в зависимости от 
значений индексов цвета (пикселей) [2]. На рис. 1 6 представлен результат сканирования композитного матери- 
ала в виде растрового изображения сечения образца правильной формы (сканирование выполнено 
В. Г. Толстиковым). Обозначены реальные изменения механических характеристик волокон и матрицы (клея) в 
структуре. Кроме того, сканирование показало возможность оценить качество композитных материалов (распо- 
ложение волокон, толщина слоя клея, внутренние дефекты и др.). 
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Рис. 1. Анализ растрового изображения сканирования материала природного происхождения 
(а) и композитного материала (6) 


Общая схема математического моделирования и алгоритм предлагаемой технологии сканирования 
представлена на рис. 2 |2]. 
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Рис. 2. Математическое и конечноэлементное моделирование реального деформируемого твердого тела 





Анализ НДС КЭ модели 


Схема отражает полный комплекс модулей: от сканирования ДТТ до получения его КЭ модели и ана- 
лиза напряженно-деформированного состояния (НДС). При необходимости внутри каждого блока алгоритма 
обеспечивается цикличность операций. 

Сканирование проводится с использованием компьютерного томографа (КТ). При работе с деформиру- 
емыми объектами он позволяет формировать пакет снимков, каждый из которых представляет собой растровое 
изображение в определенном сечении [3, 4]. Стандартные форматы растрового изображения: ]рее, Ч1сот, Бтр, 
рие и др. Независимо от используемого формата растровое изображение необходимо представить как числовую 
матрицу для дальнейшей обработки. 

Чтобы определить зависимость механических характеристик (например, в виде модуля упругости) от 
значений индексов цвета, применяется специальный весовой коэффициент [1], устанавливающий зависимость 
между средним значением индексов цвета ие, и опытным (экспериментальным) значением Ёоп. Среднее значе- 
ние модуля упругости (Еоп), полученное в эксперименте на растяжение (сжатие) стандартных образцов, ис- 
пользуется при интерпретации механических изменений по пиксельной характеристике. 

Тогда весовой коэффициент, определяющий переход от индекса цвета пикселя к модулю упругости, 
определяется выражением вида: 


Е 
пн-- ОЙ, (1) 
Пер 


Информатика, вычислительная техника и управление 


> 
©.) 
>) 


БИр://уезбиК.Чоп$ва.га 


284 


Вестник Донского государственного технического университета. 2019. Т. 19, № 3. С. 281-289. 155М 1992-5980 е155 М 1992-6006 
Тези! о{ )оп 5ие Тесйтшса Отуегуйу. 2019. Ио. 19, по. 3, рр. 281-289. 155М 1992-5980 е155М 1992-6006 





Здесь и. — среднее значение индексов цвета пакета растровых изображений сканирования: 


У льл 
В (2) 


Однако для большого объема данных сканирования КТ ДТТ, имеющих высокую степень неоднородно- 
сти структуры, при определении среднего значения индекса цвета необходимо использовать зависимость мате- 
матического ожидания. То есть, и.„ определяется как среднестатистическая величина: 


пер = Е(п)= У трит), (3) 
1=1 
где и; — значение индекса цвета на интервале Г; р;(п;) — вероятность появления п; (значения индекса цвета) на 
интервале /. 

Пикселям, каждый из которых имеет соответствующий индекс цвета, присваиваются значения модуля 
упругости материала Е’. Для этого используется весовой коэффициент Ах: 

Е; =И;КК. (4) 
Здесь Ё; — значение модуля упругости, соответствующее значению номера индекса цвета и;; Ке — весовой ко- 
эффициент модуля упругости. 

Реализация зависимости (4) для КЭ модели выполняется решением следующих задач: 

— определение зависимости механических характеристик от значений индексов цвета; 
— определение механических характеристик материала в узлах и в конечных элементах КЭ модели. 

Изменение модуля упругости определяется на основе линейной или нелинейной зависимости от значе- 
ний индексов цвета. Для материала, имеющего невысокую степень неоднородности, зависимость может быть 
построена по линейному закону (4) с использованием единственного постоянного значения весового коэффи- 
циента модуля упругости. 

Однако в природе и технике такие материалы встречаются редко. Чаще всего в структуре материала 
ДТТ можно выделить несколько областей, каждая из которых имеет отличный от других диапазон изменения 
модуля упругости. 

Для таких материалов предлагается применять нелинейный закон, позволяющий использовать в каж- 
дой области структуры ДТТ свое значение весового коэффициента. В этом случае зависимость изменения мо- 
дуля упругости ДТТ от индексов цвета представляется в виде кусочно-линейной функции (рис. 3 а). 
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Рис. 3. Кусочно-линейная (а) и непрерывная (6) функция зависимости изменения модуля упругости ДТТ 
от значений индексов цвета 

Использование кусочно-линейной функции затруднительно с точки зрения эффективности вычисли- 
тельного алгоритма, в частности, по следующим причинам: 
— чем больше количество областей структуры ДТТ, тем больше количество весовых коэффициентов и тем 
сложнее вычисление модуля упругости для узлов и конечных элементов; 
—щ переходные значения индексов цвета между областями ДТТ или объектами представляют собой не одно 
значение, а диапазоны значений индексов цвета на растровом изображении (рис. 3 6). 

Для решения этой проблемы на основе кусочно-линейной функции можно построить нелинейную зави- 


симость в виде сплайна, например кубического [5]. Такая нелинейная зависимость модуля упругости ДТТ от 
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индексов цвета является универсальной при объединении нескольких зон с различным порядком механических 
характеристик (рис. 3 6). 

Далее. Интерпретация данных сканирования связана с интерполяцией и передачей значений индексов 
цвета от пикселей к узлам КЭ сетки. В самом сложном и общем случае узлы КЭ сетки лежат между соседними 
плоскостями с растровым изображением. Они существенно различаются по форме и размерам. Кроме того, в 
них появляются дополнительные внутренние контуры и другие геометрические особенности. В этом случае для 
интерполяции используется ряд математических зависимостей, определяемых относительно прямой, наклонной 
к оси / [6]. 

В [2] представлено более подобное описание математической зависимости и блок-схемы алгоритма мо- 
дуля определения механических характеристик материала в узлах и конечных элементах КЭ модели ДТТ 
(рис. 2, блок 7). 

Для использования в КЭ модели результаты определения неоднородности механических характеристик 
материала преобразуются в файл с расширением *.рс[. Для этого задействуют две функции языка Раёгап Сот- 
тапа Гапоцазе: «табета].стеже» и «@ететргорз _стеже» [7, 8]. 

2. Исследование точности и сходимости результатов анализа КЭ моделей с учетом моделирова- 
ния неоднородности механических характеристик. В рамках данной работы определены точность и сходи- 
мость численного решения методом конечных элементов (МКЭ) с применением представленной технологии. С 
этой целью использованы расчеты НДС костной ткани человека (фрагмент бедренной кости). 

Выбор ДТТ и его анализ непринципиальны, но обусловлены важными факторами: неоднородностью 
материала кости, индивидуальностью ее геометрии, а также высоким уровнем технологии и качества сканиро- 
вания КТ в медицине. Костная ткань хорошо изучена на практике, что позволяет обоснованно судить о кор- 
ректности ее математического моделирования. 

В первом исследовании используются образцы в форме параллелепипедов, вырезанные из фрагмента 
бедренной кости человека (рис. 4 а). Результаты исследования представлены на рис. 4 6 в виде доказательства 


точности и сходимости численного решения МКЭ относительно данных натурного эксперимента [9, 10]. 
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Рис. 4. Визуализация образцов: размеры и ориентация реальных образцов (а); график сходимости нормального напряжения 
в центре образцов (6), в продольном (1), окружном (2) и радиальном (3) направлениях 


Результаты показывают, что для достижения необходимой точности численного решения МКЭ 
требуется КЭ сетка плотностью три и более конечных элемента на 1 мм [11]. Очевидно также, что свойство 
неоднородности механических характеристик материала ДТТ может быть отражено набором конечных 
элементов с собственными модулями упругости и изотропной структурой материала. 
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Рис. 5. Результат анализа НДС с применением КЭ типа НЕХ (а) и ТЕТ (6); график сходимости результатов (в) 
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Для изучения применимости двух различных типов объемных КЭ (гексаэдр и тетраэдр) проведено 
дополнительное исследование с использованием несовместных функций формы (рис. 5) [12]. Результаты 
показывают, что КЭ типа гексаэдр имеет лучшие показатели по точности и сходимости результатов, а также по 
ресурсным затратам. 

3. Моделирование реальных деформируемых твердых тел. Представленная технология применена к 
реальным объектам. С этой целью построено и проанализировано НДС КЭ моделей зубов с дефектом и 
пломбой из композитного материала. В одной модели представлен дефект в виде кариозной области (впадины) 
в верхней части эмали. Во втором варианте КЭ модели эта область заполнена пломбой из композитного 
материала. Обе модели построены с учетом неоднородности механических характеристик материала и 
индивидуальной геометрии. 

Каркасная модель зуба человека, представленная на рис.ба, построена по результатам 
сканирования [13, 14]. Изменение механических характеристик материала в объемной КЭ модели (рис. 6 6) 
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смоделировано в специальной программе [ 13, 14]. Входные данные: результаты сканирования и натурных 
испытаний. 
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Рис. 6. Построение и анализ КЭ моделей зуба человека с учетом неоднородности механических характеристик материала: 
каркасная модель (а); изменение модуля упругости зуба в объемной КЭ модели (6); расчетная схема (в); твердотельная 
геометрическая модель (г); генерация КЭ сетки (д); результат анализа НДС КЭ модели зуба с моделированием дефекта (е); 
результат анализа НДС КЭ модели зуба с моделированием композитной пломбы (ж) 


Расчетная схема представлена на рис. 6 в. На верхнем зубе моделируются дефект и композитная 
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пломба (рис. бгид). 

Результат анализа КЭ модели с дефектом (рис. бе) показывает, что кариозная часть зуба не 
воспринимает нагрузку. То есть на здоровую часть зуба приходится дополнительная нагрузка, и вероятность 
разрушения увеличивается. Результат анализа второй КЭ модели, после заполнения области кариеса пломбой 
из композитного материала (рис. 6 ж), показывает, что зуб практически полностью восстанавливает свои 
функциональные способности, так как картина напряжений идентична картине, полученной для здорового 
зуба [2]. 

Из рис. 6 ж видно, что максимумы напряжений наблюдаются в точках контакта на вершинах зубных 
бугров. Более высокая нагрузка воспринимается областью эмали и дентина. Внутренняя часть зуба менее 
нагружена. Данный результат сравнивается с расчетом [2]. При этом используется усредненное значение меха- 
нических характеристик эмали и дентина. Неоднородность структуры материала зуба не учитывается. Сравне- 
ние показывает более равномерное распределение напряжений в теле зуба при моделировании неоднородности 
распределения его механических характеристик. 

Обсуждение и заключения. Результаты исследования подтверждают точность и надежность представ- 
ленного метода КЭ моделирования на основе сканирования реальных ДТТ с учетом неоднородности механиче- 
ских характеристик материала. 

Представленный комплекс математических методов моделирования интерпретации растровых изобра- 
жений сканирования КТ позволяет проводить исследование любых сложных структур реальных ДТТ. Результа- 
ты такого моделирования используются при построении КЭ моделей с учетом неоднородности механических 
характеристик материала. 

Разработанная технология математического моделирования может быть применена для любых физиче- 
ских принципов сканирования (рентген, ультразвук, лазер и др.) и для любых типов материалов, если получен- 
ные данные представляют собой цифровое (растровое) изображение. 
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